Quench Dynamics of Topological Maximally-Entangled States 
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We investigate the quench dynamics of the one-particle entanglement spectra (OPES) for systems 
with topologically nontrivial phases. By using dimerized chains as an example, it is demonstrated 
that the evolution of OPES for the quenched bi-partite systems is governed by an effective Hamil- 
tonian which is characterized by a pseudo spin in a time-dependent pseudo magnetic field S(k,t). 
The existence and evolution of the topological maximally-entangled edge states are determined by 
the winding number of S(k,t) in the fc-space. In particular, the maximally-entangled edge states 
survive only if nontrivial Berry phases are induced by the winding of S(k,t). In the infinite time 
limit the equilibrium OPES can be determined by an effective time-independent pseudo magnetic 
field S e g-(fc). Furthermore, when maximally-entangled edge states are unstable, they are destroyed 
by quasiparticles within a characteristic timescale in proportional to the system size. 
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Topological phases that are not characterized by local 
order parameters have been subjects of key interests in 
condensed matter physics due to anomalous properties 
associated with these phases. Recent discovery of time- 
reversal invariant topological insulators 0, 0] has further 
triggered intense investigations on the characterization 
of topological phases. One of the important features for 
topological phases is the possibility of creating nonlocal 
properties associated with the topology, which is often 
realized as the entanglement between the system and its 
environment. In particular, the entanglement spectrum 
(ES) i.e. the eigenvalues of the reduced density matrix 
of the system, provides an idea tool to characterize the 
topological phase @, 3 . In general, not only the ES can 
be used to distinguish different classification of topolog- 
ical phases but also it can detect the existence of edge 
modes at zero energy Q . The existence of edge states 
at zero energy reflects the non-trivial topology of the un- 
derlying quantum state. From the point of view for ma- 
nipulating quantum information, these edge modes rep- 
resent the maximally entangled states of the system and 
its environment. Therefore, they are the best candidate 
for qubits @. 

In order for the topological maximally-entangled state 
being a viable candidate for qubits, it is necessary to ex- 
amine if they could survive under quantum information 
processing. Since typical quantum manipulations involve 
rapid change of the coupling to the environments, it is 
therefore important to examine the quench dynamics of 
maximally-entangled states. Recent investigations indi- 
cate that the thermalization of integrablc systems due 
to quench depends strongly on the initial conditions 0- 
[Toj | . These studies, however, are confined to the bulk 
properties. There are only few papers concerning quench 
dynamics of topological edge states [ll| . In the presence 



of quench dynamics could be entirely different. It is thus 
important to examine the quench dynamics of topologi- 
cal maximally-entangled states. In this paper, by taking 
dimerized chain as an example, we investigate how edge 
states affect the thermalization and the quench dynam- 
ics of one-particle entanglement spectra (OPES) defined 
below. In particular, we show that the existence and 
evolution of topological maximally-entangled edge states 
arc determined by the winding number of a pseudo mag- 
netic field S(k,t). The maximally-entangled states sur- 
vive only if nontrivial Berry phases are induced by the 
winding of S(k,t). 

Consider the ground state I'J'ab} of a bipartite total 
system AB that consists of the system A and the environ- 
ment B. The reduced density matrix of the system A is 
PA = Tr B |4'ab}(^ , ab|- The entanglement entropy (EE), 
defined as 5a = — TrpA log 2 Pa, has been widely used to 
measure the bipartite entan glem ent between the system 
A and the environment B 12]. It is known that the 
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of topological edge states, the system can be maximally 
entangled with the environment and the thermalization 



scaling law of EE provides a way to distinguish different 
quantum phases [131 ] . Furthermore, the property that EE 
diverges at the critical points provides a useful tool to ex- 
amine the quantum criticality [141 ] . In addition to global 
properties associated with EE, it is useful to explore de- 
tailed microscopic quantum phenomena using OPES, de- 
fined as the set of A m 's with pa — ® m J? , ^ -> 

U 1 — A. 

and OPES arc related through the relation Sa = 
where S m = -A m log 2 A m - (1 - A m ) log 2 (1 - A m ). The 
OPES has been used to investigate disorder lines flij ]. 
Berry phase @, Q and zero-energy edge states @, 0] ■ It 
is clear that the eigenvalue A m = 1/2 corresponds to the 
situation when the system A and the environment B are 
maximally entangled so that S m = 1. Since A m = 1/2 is 
a mid-gap state that often results from the zero-energy 
edge state between A and B, the maximally- entangle 
state is often topologically protected. 
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To investigate quench dynamics of the maximally- 
entangled state, we consider a one-dimensional (ID) 
dimerized chain characterized by the Hamiltonian: 



(1) 



where i is the site index and <5 £ [—1,1]. The ground 
state undergoes a phase transition from a topologically 
trivial phase (6 < 0) to a topologically nontrivial phase 
(S > 0) as S is varied across the phase boundary S c = 0. 
Let the region i > represents the system A while i < 
represents the environment B, the topological edge mode 
appears if the bond [i, i + 1] = [0, 1] is strong, i.e., —(1 + 
6) < —(1 — 6). Otherwise the edge mode does not exist 
1 1 51 ] . In this paper we study the quench dynamics of the 
topological maximally-entangled states for a dimerized 
chain by suddenly quenching the parameter 5. 

The occurrence of maximally-entangled edge states 
has its topological origin. By defining a spinor c Xi — 
{cxn , Cx i2 ) T where xn = 2i — 1 and x^ = 2i and per- 
forming a Fourier transformation, the Hamiltonian (TTJ) 
can be castcd into the form 



H s = - c I[ R i( fc ) -o-]cfc. 



(2) 
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where a = (cr x ,(Ty,cr z ) arc Pauli matrices and 

R 5 (k) = ((1 — 6) + (1 + 6) cos k, (1 + 8) sin k, 0) (3) 

is a pseudo magnetic field with the magnitude Rs(k) = 
2(cos 2 fc/2 + S 2 sm 2 fc/2) 1 / 2 . The system fulfills chiral 
symmetry because R,s(A;) lies on a plane. For 8 £ (0, 1] 
the loop £ of Rs(k) encloses the origin O as k runs 
through the Brillouin zone. Consequently, the Berry 
phase (or Zak's phase) Xi defined as a line integral of the 
curvature of the filled band, is it. Due to the fact that £ 
can be continuously deformed into a unit circle without 
crossing the origin, topological argument ensures that the 
original Hamiltonian corresponding to I contains at least 
one pair of one zero-energy edge states as a consequence 
of chiral symmetry On the other hand, if 8 £ [—1,0), 
one obtains trivial Berry phases and no topological edge 
state occurs (sec the upper panels of Fig[T]). 

We consider a total system of infinite size, and parti- 
tion it into the system A with Xi = 1, • • • , L/2 and its 
environment B (the remaining part). The reduced den- 
sity matrix of the system pA for the ground state can 
be determined by eigenvalues X m of the block Green's 
function matrix (GFM) O G s a p(xi,Xj) = Trfpct.^.J, 
where Xi a ,Xj@ belong to the system A, a,/3 = 1,2, and 
p is the ground state density matrix of the total system. 
Hence the GFM can be considered as an effective Hamil- 
tonian that determines the OPES. In the Fourier fc-spacc, 
G a/3( fc ) is g iven b y 
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FIG. 1: (color online) Upper two panels: Representative loops 
I of R^(fc) for phase I and II with 8 — 0.5 and 5 — —0.5 
respectively. In phase I (<5 6 (0, 1]) the Berry phase is n 
due to circling of R^(fc) around the origin, while in phase II 
(8 6 [—1,0)) the Berry phase is trivial. Lower four panels: 
evolutions of two A m (t) that are closest to 1/2 for different 
quench processes: (a) phase I (8 = 0.5) to I (8 — 0.2) (b) 
phase I (8 = 0.5) to II (8 = -0.5) (c) phase II (8 = -0.5) to 
I [8 = 0.5) (d) phase II (8 = -0.5) to II (8 = -0.2). Insets: 



1 + R«(fe) a- 



a/3 



(4) 



loop I of S e g- 



where Ra(fc) = R^(fe)/ 'Rs(k). It takes almost the same 
form as the Hamiltonian (J2|) except for a constant and 
a positive normalization factor R$, which leads to the 
conclusion that they share the same topology. If R^(fc) 
encloses the origin in the parameter space with Berry 
phase equal to n, a pair of zero energy states appear for 
the Hamiltonian @, while for the GFM, we obtain a pair 
of maximally-entangled states with A rn = 1/2. Notice 
that the maximally-entangled states are also edge states, 
which can be seen in Fig [31 

Consider now a sudden quench at t = by changing 
the parameter of the dimerized chain in Eq. (TTJ) from 
5 at t < to 5' at t > 0. Denote the phase whose 
S £ (0.1] with maximally-entangled states by phase I, and 
the phase whose 6 £ [—1, 0) by phase II. We perform four 
possible quench processes, including I to I, I to II, II to I, 
and II to II. The OPES A TO (t) at time t can be obtained 
by diagonalizing a time-dependent GFM which is de- 
fined as G a 3(xi,Xi,t) = Tr(/9e iH «'*4.. Q C x ^e" li L'*) where 
Xi,xj £ A [18(. In Fig[T] (a)-(d), we show the time evolu- 
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tion of OPES for states whose eigenvalues are closest to 
1/2. In cases of I to I and I to II, the maximally-entangled 
states persist for a while and then split into two differ- 
ent evolutions. Only for the case of I to I, however, the 
two splitting eigenstates evolve back to two maximally- 
entangled states at infinite time. In contrast, for the case 
of I to II, two splitting eigenstates closest to 1/2 remain 
splitting forever. On the other hand, if one starts with 
initial states without maximally-entangled states such as 
the cases of II to I and II to II, maximally-entangled 
states can not be created at later time. 

It is instructive to define a time-dependent pscudo 
magnetic field S(k,t) from the time-dependent GFM 
G(k,t) in the Fourier space through the relation 



G a/3 (M) = l/2 [1 + S(MWL 



/3 • 



(5) 



S(k,t) can be further written as a summation of three 
different contributions 

Si(M) = cos(2Rs>t)R s (k), (6) 
S 2 (M) = sm{2R s d)R s (k) xR 5 ,(k), (7) 
S 3 (M) = (1 - cos{2R 5 d))R s (k) • R*,(fc)Rj/(fc). (8) 

The thermalization of EE observed in Fig[T]can be under- 
stood by considering the infinite-time limit of the time- 
dependent pseudo magnetic field ([5]), ©, and (J5J). Since 
the sinusoidal parts dephase out, the long-time behav- 
ior of the time-dependent GFM, G a p{k,t = oo), is solely 
determined by the effective pseudo magnetic field 



S cff (fc) 



R s (k)-R s ,(k))R s ,(k). 



(9) 



The existence of maximally-entangled states at infinite 
time is thus determined by the topology of S c g\ If 
S c g- encircles origin, the Berry phase is tt, maximally- 
entangled states appear at t = oo; otherwise there is 
no maximally-entangled state at later time. It is clear 
that both the final Rs' and the initial Rs determine S e Q- 
It implies that the long time behavior carries the mem- 
ory of the initial state, which is due to nonergodicity 
of integrable systems The existence of maximally- 

entangled states (edge states) at infinite time requires 
that both the initial and final Hamiltonians possess non- 
trivial Berry phase. In the inset of Fig HJa)-(d), S g fj 
are plotted for different quench processes with different 
initial and final Hamiltonians. It is clear that only for 
the quench from I to I, S e g- encircles origin, while in the 
other quenches S e -jj either not encircles the origin (II to 
II) or passes through the origin (I to II or II to I). 

The reason is that if the initial state is in the same 
phase as the final Hamiltonian (I to I or II to II), 
Rs{k) ■ Rs'(k) is always positive, therefore the topology 
of S c g- is the same as Rs'(k). On the other hand, if the 
system is quenched into a phase that is distinctly from 
the initial state, Rg(k = 0) • Rg'(k = 0) is positive while 
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FIG. 2: (color online) Evolutions of two X m (t) that are closest 
to 1/2 for different L under the quench of I (8 = 0.5) to II 
(S = —0.5). Lower inset: X m (t) obtained from exact diago- 
nalization of G ctl ^(xi,Xj,t — 1000) (blue circle) and S e g- (red 
dots). Upper left panel: fitting of T* , the persistence time of 
the maximum entangled states, using T*(L) = a(L)/(l — \5'\) 
where T* is the persistence time of the maximum entangled 
states. Upper right panel: fitting of a(L) versus L. We find 
that a « 0.232i. 



Rg(k = tt) ■ R$> (fc = tt) is negative, then there must exist 
one point fco that S e fj(fco) ^ R-s(ko) -Rg' (ho) = 0. Hence, 
ko destroys the topology of the infinite-time Green's func- 
tion matrix. There exist no maximally-entangled states 
whenever the system is quenched across the topological 
phase boundary. This explains why only the quench from 
phase I to phase I creates maximally-entangled states at 
infinite time. 

For all cases shown in Fig [TJ the OPES A m (t) of the 
dimcrized chain fluctuate before it reaches the equilib- 
rium. The intermediate regions can be explained by the 
appearance of S 2 in the time-dependent GFM. Since S 2 
is proportional to R$(k) x R$r(k) and is perpendicular 
to Si and S3, the system is agitated by quasiparticles in- 
duced by S 2 and S3 until the time-dependent sinusoidal 
functions dephase out and then the system reaches its 
equilibrate state with the effective S e g- ([9]). Clearly, the 
thermalization depends on sizes of the system A. To 
check the size dependence. Fig [2] shows several quench 
processes from phase I to phase II with different sys- 
tem length. In the lower panel, we compare eigenval- 
ues at t = 1000 (red dots) evolved from the maximally- 
entangled state with eigenvalues obtained by diagonal- 
ization of G Q< 5(k, 00) (blue circles). The dependence on 
the system size L shows excellent agreement, indicating 
the validity of S c g-. 

We also find that if the initially the system is in phase 
I, the maximally-entangled states can persist for a long 
time T* before they are destroyed. This feature is strik- 
ing because the maximally-entangled states are not de- 
stroyed in the beginning. Furthermore T* is independent 
of the initial conditions but depends on the final Hamil- 
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FIG. 3: The probability sum of the two eigenstates whose 
eigenvalues are closest to 1/2. (a) Quench process I to I. The 
edge modes persist from the beginning to the end. (b) Quench 
process I to II. The edge modes disappear and evolve to the 
bulk state. 



tonian and the system length L. This is because the 
maximally-entangled states reside on the edge Q and 
they will disappear only when the quasiparticles created 
from the bulk reach the edges. Therefore, we expect 
T* re L/(2v m a,x), where fmax is the maximum veloc- 
ity of the quasiparticles. For the energy dispersion on 
obtains Umax = 2(1 — \5'\) hence 



4(1-|«J'|)- 



(10) 



In the upper left inset of Fig [2J we show the fitting of 
T*(L) to the function a(L)/(l — 1<5'|) using various system 
size L. We then fit a{L) = b x L as shown in the upper 
right inset of Fig [2] to find b re 0.232. Combining these 
two fitting we find that T* re 0.232L/(1 - \S'\) which is 
very closed to our approximation Eq. (fit))) . Our results 
imply that edge modes remains to the infinite time if L is 
infinite. Therefore, the edge modes serve as good qubits 
since they have the maximal entanglements and cannot 
be destroyed easily if the system is large enough. 

The disappearance and recreation for different quench 
process starting from the phase I can be further explored 
by examining the probability sum of the two eigenstates 
with eigenvalues being closest to 1/2. Fig EJa) shows 
the case of I to I. It is seen that only small amounts of 
quasiparticles are excited. As a result, the edge modes 
remain from the beginning to the end. On the other 
hand, for the quench process I to II as shown in Fig[3Jb), 
quasiparticles strongly modify two edge modes until they 
disappear and the whole system is bulk-like without any 
edge mode. 

In summary, using dimerized chains as an example, the 
quench dynamics of the maximally-entangled states is in- 
vestigated by diagonalizing the time-dependent Green's 



function matrix. We find that the existence of the 
maximally-entangled states after sudden quench is de- 
termined by an effective pseudo magnetic field (S e fj), 
which depends on both the initial and final Hamiltonians. 
The topological properties at infinite time are thus de- 
termined by the initial states and the Hamiltonian after 
the quench. When the maximally-entangled states are 
unstable, they are destroyed by quasiparticles that move 
from the bulk to the edges with a characteristic timescale 
proportional to the system size. 

We thank Dr. Sungkit Yip for useful discussions. This 
work was supported by the National Science Council of 
Taiwan. 
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